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The Abundance Pattern of O, Ne, Mg, and Fe in the Interstellar 
Medium of SO Galaxy NGC 4382 Observed with Suzaku 
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Abstract 

We derived O, Ne, and Mg abundances in the interstellar medium (ISM) of a relatively isolated SO 
galaxy, NGC 4382, observed with the Suzaku XIS instruments and compared the O/Ne/Mg/Fe abundance 
pattern to those of the ISM in elliptical galaxies. The derived temperature and Fe abundance in the ISM 
are about 0.3 keV and 0.6-2.9 solar, respectively. The abundance ratios are derived with a better accuracy 
than the abundances themselves: O/Fe, Ne/Fe, and Mg/Fe ratios are 0.3, 0.7, and 0.6, respectively, in 
solar units. The O/Fe ratio is smaller than that of the ISM in elliptical galaxies, NGC 720, NGC 1399, 
NGC 1404, and NGC 4636, observed with Suzaku. Since O, Ne, and Mg are predominantly synthesized 
by supernovae (SNe) of type II, the observed abundance pattern indicates that the contribution of SN la 
products is higher in the SO galaxy than in the elliptical galaxies. Since the hot ISM in early-type galaxies 
is an accumulation of stellar mass and SN la products, the low O/Fe ratio in the ISM of NGC 4382 reflects 
a higher rate of present SNe la, or stars containing more SN la products than those in elliptical galaxies. 
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1. Introduction 

Early- type galaxies have a hot, X-ray emitting inter- 
stellar medium (ISM), which is considered to be gravita- 
tionally confined (e.g., Forman et al. 1985; Mathews <fc 
Brighcnti 2003). X-ray observations of the metal abun- 
dances in the ISM of early-type galaxies provide a key to 
understanding the history of star formation and the evo- 
lution of galaxies, because the metals in the ISM come 
from type la supernovae (SNe la) and stellar mass loss. 

A lot of observational evidence suggests that a signif- 
icant fraction of present early-type galaxies have trans- 
formed from late-type galaxies (e.g., Dressier et al. 1997; 
Fasano et al. 2000; Treu et al. 2003; Postman et al. 
2005; Smith et al. 2005). From z-0.5, the fraction of 
SO galaxies in clusters of galaxies has increased, whereas, 
that of spirals has decreased, with no evolution in the frac- 
tion of ellipticals (Kodama et al. 2004; Desai et al. 2007). 
Poggianti et al. (2009) found that these changes appear 
more strongly in less massive clusters with lower veloc- 
ity dispersion. These results suggest that spiral galaxies 
changed into SO galaxies at z<0.5, falling into clusters of 
galaxies. The mctallicity of stars in galaxies reflects in 
the star formation history, therefore, it is an important 
parameter for understanding the evolution of galaxies. 

Optical observations indicate that Mg/Fe and a/Fe ra- 
tios of stars are super-solar in the cores of bright early- type 
galaxies and increases with galactic mass (e.g., Faber et 
al. 1992; Worthey et al. 1992; Nelan et al. 2005; Thomas 
et al. 2005; Graves et al. 2007). This overabundance of 
Mg relative to Fe is the key indicator that galaxy for- 
mation occurred before a substantial number of SNe la 
could explode and contribute to lowering the these ratios 



(e.g., Bcrnardi et al. 2003; Nelan et al. 2005; Smith et al. 
2006; Graves et al. 2007; Pipino et al. 2009). However, 
absorption-line indices that account for abundance ratios 
also depend on the age distribution of stars. Optical spec- 
troscopy is limited within the very center of galaxies. 

Using X-ray observations, we can directly determine the 
metal abundances of the ISM, and constrain the stellar 
mctallicity of the entire galaxy. The atomic data for lines 
at X-ray wavelengths are simpler than for those in op- 
tical spectra, and the structure of the hot ISM is also 
much simpler than stellar population data. Therefore, 
we can estimate the temperature and metallicity of the 
hot ISM through X-ray spectra with small systematic 
uncertainties. XMM-Newton EPIC and RGS provided 
the means to measure O and Mg abundances in some 
systems, but reliable results have been obtained only 
for several central galaxies in groups and clusters (e.g., 
Xu et al. 2002; Matsushita et al. 2003; Tamura et al. 
2003; Matsushita et al. 2007a). The ISM in such re- 
gions might be polluted by the gas of clusters or groups 
(e.g., Matsushita 2001; Matsushita et al. 2002; Nagino & 
Matsushita 2009). 

With the Suzaku X-ray satellite, O, Ne, and Mg abun- 
dances of four elliptical galaxies, NGC 720 (Tawara et 
al. 2008), NGC 1399 and NGC 1404 (Matsushita et al. 
2007b), and NGC 4636 (Hayashi et al. 2009), have been 
measured. The XIS detector onboard Suzaku (Mitsuda et 
al. 2007) can constrain O, Ne, and Mg abundances well 
because its energy resolution is better, and its background 
is lower than any previous X-ray CCD detector (Koyama 
et al. 2007). According to the new solar abundance table 
of Loddcrs (2003), the abundance ratios of these elliptical 
galaxies are close to the solar values. Assuming the SN II 
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abundance pattern of Iwamoto et al. (1999), about 80% 
of Fe is synthesized by SNe la (Matsushita et al. 2007b). 

The abundance pattern of the hot ISM in SO galaxies 
observed with Suzaku has not been reported. In this pa- 
per, we present the temperature and abundances of O, 
Ne, Mg, and Fe in the ISM within 4 r e of the SO galaxy 
NGC 4382, observed with Suzaku. Here, r e is the effec- 
tive radius of the galaxy. For NGC 4382, r e corresponds 
to 0.91 arcmin (de Vaucouleurs et al. 1991). The distance 
to NGC 4382 is 16.8 Mpc (Tally 1988), and its redshift, 
z=0.002432, is taken from the NASA/IPAC extragalac- 
tic database (NED). NGC 4382 is located in the outskirts 
of the Virgo cluster, 1.7 Mpc from cD galaxy M87. The 
central stellar velocity dispersion of the galaxy, 179 km 
s -1 and its [a/Fe] value, 0.12 ±0.06 (McDermid et al. 
2006), are much smaller than those of the four elliptical 
galaxies, observed with Suzaku. These values were de- 
rived only within two arcsec of the center of NGC 4382, 
and we do not know the metallicity of stars outside of the 
central region. So far, NGC 4382 has been observed with 
ROSAT (Fabbiano et al. 1994), ASCA (Kim et al. 1996), 
Chandra (Sivakoff et al. 2003; Athey 2007) and XMM- 
Newton (Sansom et al. 2006; Nagino & Matsushita 2009). 
These X-ray observations suggest that the temperature of 
the ISM in NGC 4382 is about 0.3-0.4 keV, which is much 
lower than those of the four early-type galaxies observed 
with Suzaku. Since no the intracluster medium (ICM) 
around NGC 4382 has been detected, this galaxy is suit- 
able for investigating heavy elements in the ISM of the 
galaxy itself. 

Throughout this paper, we use a Hubble constant of 
Hq = 70 km s _1 Mpc -1 . We adopt the new solar abun- 
dance table of Lodders (2003). The abundances of O and 
Fe have increased by about 70% and 60%, respectively. 
Unless otherwise specified, errors arc quoted at 90% con- 
fidence. 

2. Observation and data reduction 

Suzaku observed NGC 4382 in June 2008 with an ex- 
posure of 99 kscc in the XIS nominal position, pointing 
at the center of NGC 4382. Figure 1 shows a 0.3-2.0 keV 
image of NGC 4382. We analyzed only the XIS data in 
this paper, although Suzaku observations were performed 
with both XIS and HXD. The XIS instrument consists of 
three sets of X-ray CCDs (XIS0, 1, and 3). XIS1 is a 
back-illuminated (BI) sensor, while XIS0 and 3 are front- 
illuminated (FI) ones. In this observation, these instru- 
ments were operated in the normal clocking mode (8 s 
exposure per frame), with the standard 5 x 5 or 3 x 3 edit- 
ing mode. 

Our analysis was performed using XSPEC version 
11.3.2ag and HEAsoft version 6.6.2. We processed the XIS 
data with xispi, using the calibration files distributed on 
04-03-2009. Then, we cleaned the XIS data by assuming 
thresholds of > 5° on the Earth elevation angle, and > 20° 
on the day-earth elevation angle. The total exposure time 
for both FI and BI sensors amounts to ~99 ksec after stan- 
dard data screening. Event screening with cut-off rigidity 




Fig. 1. Raw XIS image of NGC 4382 within an energy range 
of 0.3—2.0 keV. XIS0, 1, and 3 images were added in sky co- 
ordinates. Non- X-ray background (NXB) and cosmic X-ray 
background (CXB) were not subtracted. Image was smoothed 
with a Gaussian of a = 0.56 arcmin. White circles represent 
regions within 4 r e and 5 r e , respectively. Point source indi- 
cated by magenta circle is excluded from spectra. 

(COR) was not performed on our data. 

XIS response files were calculated using the xisrmf gen 
response matrix file (RMF) generator (Ishisaki et al. 
2007), version 02-28-2009. Ancillary response files (ARF) 
were calculated using xissimarf gen, version 01-08-2009, 
assuming flat-sky emission. We note that ARF files as- 
suming a point-source, a /3-model profile, or flat-sky are 
almost the same within the energy range of 0.4 to 5 keV, 
except for normalization. Slight degradation of the energy 
resolution was considered in the RMF, and a decrease in 
the low-energy transmission of the XIS optical blocking 
filter (OBF) was included in the ARFs. To subtract the 
contribution of non- X-ray background (NXB), we used 
xisnxbgen, which generates the NXB spectrum from a 
dark Earth database (Tawa et al. 2008). 

3. Spectral analysis 

As shown in Figure 1 , we extracted the spectra within a 
4 r e circle centered on the optical center of NGC 4382 from 
NED. We also used the region outside 5 r e to estimate 
the background. Northwest of NGC 4382, a point-like X- 
ray source was detected at 0.3-2.0keV. We excluded this 
source from each spectrum. 

We used the spectra from the BI (XIS1) and FI (XIS0 
and 3) sensors for the energy ranges of 0.4-5.0 keV and 
0.6-5.0 keV, respectively, since subtracting background 
lines is difficult above 5.0 keV, and the ISM does not 
emit photons above 5.0 keV. The energy range around 
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the Si K-edge (1.82-1.84 keV) was ignored as a result 
of a problem in the response matrix. The remaining to- 
tal photon counts with XISO, 1, and 3 are 3490, 6604, 
and 3617 counts, respectively. Although two calibration 
sources emitting strong lines above 5 keV are located at 
the corners of each XIS, spectral analysis is not affected 
below 5 keV. Therefore, we included these regions, and 
analyzed the spectra below 5 keV. The spectra of XIS de- 
tectors (XISO, 1, and 3) were fitted simultaneously with a 
model consisting of the ISM emission, unresolved discrete 
sources, cosmic X-ray background (CXB), and Galactic 
components. Here, the spectra of the NXB were sub- 
tracted. 

3.1. The ISM component 

We adopted a single-temperature vAPEC model (Smith 
et al. 2001) as the ISM thin thermal emission component 
(hereafter, the IT model). The ISM in early-type galaxies 
originates mainly from stellar mass loss, and most of the 
metals in these stars are the products of SNe. The metal 
abundances of He, C, N, and Al were fixed at solar values. 
We organized the abundances of other heavy elements into 
five groups: O; Ne; Mg; the Si group (Si, S, Ar, and Ca); 
and the Fe group (Fe and Ni), and allowed them to vary. 
This is because significant detections of the ionized O, 
Ne, Mg, Si, and Fe lines are expected. In this study, we 
assume that abundances of S, Ar, and Ca are the same as 
that of Si, taking into account the metal synthesis process 
in SNe. These lines were not significantly observed, and 
hence we could not determine the abundance of each line. 
In our Galaxy, Ni/Fe ratios synthesized from both SN la 
and SN II are solar values (Edvardsson et al. 1993). 

This component was subjected to a common absorption 
with a fixed column density (iVn) at the Galactic value of 
Dickey & Lockman (1990), 2.5 x 10 20 cnT 2 . 

3.2. Unresolved discrete sources 

For unresolved point sources, we added a power-law 
model, in which the photon index was fixed at 1.6 (Irwin et 
al. 2003). The power-law model described the total spec- 
tra of discrete sources in early-type galaxies well, and was 
successfully used for this role in previous studies (Blanton 
et al. 2001; Randall et al. 2004). This component, as well 
as the ISM component, was also subjected to a common 
absorption with fixed Nn at the Galactic value of Dickey 
& Lockman (1990). 

3.3. CXB and Galactic emission 

For the CXB component, we adopted a power- law 
model with a photon index L = 1.4 (Kushino ct al. 2002). 
Galactic emission arises mainly from the local hot bubble 
(LHB) and the Milky Way halo (MWH). We used a two- 
temperature APEC model for these components and fixed 
the temperatures of the LHB and MWH emission at 0.1 
and 0.3 keV, respectively, following previous studies (e.g., 
Komiyama ct al. 2009; Konami et al. 2009; Hayashi ct al. 
2009; Yamasaki et al. 2009). CXB and MWH components 
were subjected to absorption at the Galactic value, while 
the LHB was not. 
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Fig. 2. Spectra outside of 5 r e observed with XISO (black), 
XIS1 (red), and XIS3 (blue). These spectra are fitted with 
a sum of the CXB and Galactic models. Lower panel shows 
fit residuals. Here, for simplicity, only model components for 
the XIS1 spectra are shown. Green and gray lines correspond 
to CXB and Galactic emission components, respectively. 

As shown in Figure 1, we fitted the spectra with the 
NXB subtracted outside of 5 r e from the center of NGC 
4382 to estimate the contribution from the CXB and 
Galactic emission. In this study, we assumed that there 
is no emission from the galaxy beyond 5 r e . The normal- 
izations of CXB and Galactic components were allowed to 
vary. Figure 2 shows the best-fit model and spectra for 
the CXB and Galactic components. The spectra are well 
described by the model with a reduced-% 2 value of 1.23: 
there were 419 degrees of freedom. 

3.4- Spectral fitting within 4 T e 

Next, we fitted the spectra within 4 r e and outside of 
5 r e of XISO, 1, and 3 simultaneously with the IT model 
for the ISM, two power-law models for unresolved sources 
and CXB, and two APEC models for Galactic emission. 
Here, we assumed that the normalizations of the ISM and 
unresolved source components were fixed at beyond 5 r e . 
We also assumed that the CXB and Galactic emission 
within 4 r e and beyond 5 r e have the same surface bright- 
ness and spectral parameters. Figure 3 shows the XIS 
spectra fitted in this way. The derived x 2 , ISM temper- 
ature, and abundances are summarized in Table 1. This 
fit provides a rather small reduced-x 2 of ~ 1.15, although 
residual structures appear around 0.8 keV. 

4. Results 

4-1. ISM temperature and abundances 

Table 1 summarizes the best-fit parameters of temper- 
ature and abundances. The ISM temperature within 4 r e 
is about 0.3 keV. This temperature is lower than those 
of early-type galaxies, NGC 720, NGC 1399, NGC 1404, 
and NGC 4636, which have been observed with Suzaku 
(Matsushita et al. 2007b; Tawara et al. 2008; Hayashi et 
al. 2009). This temperature is consistent with those de- 
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Fig. 3. Spectra within 4 r e observed with XISO (black), XIS1 (red), and XIS3 (blue). These spectra are fitted with lT+power-law 
(left), and lT+bremss (right) models. Lower panels show fit residuals. For simplicity, only model components for XIS1 spectra are 
shown. Magenta and orange lines represent ISM components and unresolved sources, respectively. Green and gray lines correspond 
to CXB and Galactic emission components, respectively. 



Table 


1. Results of spectral fittings 


within 4 r e 


of NGC 4382 with IT models 




model 


kT 

(keV) (solar) 


Ne 
(solar) 


Mg 

(solar) 


Fe 
(solar) 


x'Vd-O.f. 


lT+power-law 
lT+brcmss 


n 011+0.1)24 n oo+U.BV 
U.01J._q 01 g U.OO_g 16 

0.30818:88 0.30185 


n 7n+ l ab 

u - ' u -0.31 

64+ 95 


63 +iyy 
u.oo_ 28 


93+ 134 


759 / 658 
754 / 658 



model 


O/Fe 
(solar) 


Ne/Fe 
(solar) 


Mg/Fe 
(solar) 


lT+power-law 
lT+bremss 


0.33+8ii 
0.32l o ;i 8 o 


"•<;»• 


0.6018:2? 



rived from previous observations with ROSAT (Fabbiano 
et al. 1994), Chandra (Sivakoff et al. 2003; Athey 2007), 
and XMM-Newton (Nagino & Matsushita 2009). On the 
other hand, Sansom et al. (2006) detected a higher ISM 
temperature of about 0.4 keV with XMM-Newton EPIC. 

The abundances are less constrained, and Fe abundance 
within 4 r e is ~ 0.6-2.9 solar with large errors. This value 
is consistent with those of the early-type galaxies, NGC 
720, NGC 1399, NGC 1404, and NGC 4636. The other 
abundances, O, Ne, and Mg, are mostly consistent with 
the Fe abundance within large error bars, but tend to have 
lower values than that of Fe. 

4-2. Abundance ratios 

Figure 4 shows the confidence contours of O, Ne, and 
Mg abundances against Fe abundance. Fe in the ISM orig- 
inates from type la and type II SNe, whereas, O, Ne, and 
Mg come from type II SNe (Iwamoto et al. 1999; Nomoto 
et al. 2006). Thus, 0-, Ne-, and Mg-to-Fe abundance 
ratios are particularly important for constraining the con- 
tributions from type la and II SNe. These contours show 
elongated shapes, especially for the Ne and O abundances 
vs. Fe abundance, and are along the lines representing 
constant O/Fe, Ne/Fe, and Mg/Fe ratios. Thus, we can 
constrain the metal-to-Fe abundance ratios well, although 



the individual error bars for each metal abundance are 
large. 

O/Fc, Ne/Fe and Mg/Fe ratios derived from these 
confidence contours arc also summarized in Table 1. 
Abundance ratios of O/Fe, Ne/Fe, and Mg/Fe in solar 
units are 0.3318:!?, 0.69l8;if, and 0.62±°:^, respectively. 
O abundance is especially small relative to the Fe abun- 
dance. 

4-3. Uncertainties in spectral fit 

In this section, we examine the uncertainties in the 
abundance determination because of an uncertainty in the 
spectral modeling. 

To investigate the systematic effect of the model for 
unresolved sources, described in Section 3.2, we also fit- 
ted the components with a bremsstrahlung model (e.g., 
Matsushita et al. 1994; Irwin et al. 2000) and compared 
the results with those using the power-law model. Here, 
we fixed the temperature of the bremsstrahlung model at 
7.0 keV (Matsushita et al. 1994). The best-fit spectra 
are shown in Figure 3. As shown in Table 1, the derived 
results are almost the same. 

We also tried a spectral fitting with a two-temperature 
model for the ISM (hereafter, the 2T model). However, 
the values of reduced-;^ 2 did not improve, and the same 
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Fig. 4. Confidence contours for O, Ne, and Mg abundances vs. Fe abundance. Crosses denote the best-fit, and black-solid, 
gray-solid, and gray-dotted contours represent 68%, 90%, and 99% confidence levels, respectively. Dashed lines correspond to 
constant abundance ratios. 



residuals remained around 0.8 keV, at the Fe-L line ener- 
gies. The temperature of the first component was almost 
the same as that in the IT model, and the normalization of 
the second component was almost zero. Thus, the derived 
abundances were unchanged. We then fitted the spectra 
with a five-temperature model for the ISM (hereafter the 
5T model); we fixed the ISM temperatures at (i) 0.2, 0.4, 
0.6, 0.8, and 1.0 keV and (ii) 0.1, 0.2, 0.3, 0.4, and 0.5 
keV. 

Figure 5 shows spectra fitted with the 5T(i) and 5T(ii) 
models. The spectral fitting with the 5T(i) model pro- 
vided a worse \ 2 value than that with the IT model (Table 
2). Only the normalizations of ISM components with tem- 
peratures of 0.2 and 0.4 keV are >0, and the normaliza- 
tions of the other ISM components fall to 0. This suggests 
that temperature components above 0.6 keV do not exist 
in the ISM of NGC 4382. On the other hand, the spec- 
tral fitting with the 5T(ii) model provides the same x 2 
as the IT model (Table 2). In this fitting, only the nor- 
malizations of temperature components at 0.3, 0.4, and 
0.5 keV arc > (cf. Figure 5). The contribution of the 
0.3 keV component is the largest of the three tempera- 
ture components at about 97%. The derived abundances 



and abundance ratios did not differ from the results of IT 
spectral fitting. These 5T models even gave similar resid- 
ual structures around 0.8 keV. These discrepancies in the 
Fc-L energy range are also seen in Suzaku observations of 
NGC 720, NGC 1404, and NGC 4636, whose ISM temper- 
atures are also ~ 0.6 keV (Matsushita et al. 2007b; Tawara 
et al. 2008; Hayashi et al. 2009). Therefore, these residual 
structures are likely to be related to poorly modeled Fe-L 
lines. 

To investigate the effect of the temperatures of the 
Galactic components on the ISM abundances of NGC 
4382, we allowed the temperature of the MWH compo- 
nent to vary and re-fitted the spectra. We got almost the 
same ISM temperature and abundances. 

4- 4- Radial profiles 

To study radial discrimination, we investigated the 
spectra extracted from the different radial regions, R < 
2r e , 2r e < R < 4r e , and 3r e < R < 6r e . Here, R is the 
projected radius from the galaxy's center. In Table 3, we 
summarized the fitting results. A low signal-to-noise ra- 
tio does not constrain the abundances. We can derive the 
abundance ratios against Fe only, although each value has 
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Fig. 5. Spectra within 4 r e observed with XISO (black), XIS1 (red), and XIS3 (blue). These spectra are fitted with 5T(i)+power-law 
(left), and 5T(ii)+power-law (right) models. Lower panels show fit residuals. For simplicity, only model components for XIS1 spectra 
are shown. Colors have the same meanings as in Figure 3 



Table 2. Results of spectral fittings within 4 r e of NGC 4382 with 5T models 



model 




kT 
(keV) 


O 

(solar) 


Ne 
(solar) 


Mg 
(solar) 


Fe 
(solar) 


xVd.o.f. 


5T(i)+power-law 
5T (ii) +power-law 


0.2 
0.1 


0.4, 0.6, 0.8, 1.0* 
0.2, 0.3, 0.4, 0.5* 


U - OM -0.15 
n o 9 +0.70 
U - Jz -0.12 


0.68™ 
0.74+° ; 8 25 


f| K7+t.2b 
u - ' -0.29 

fiQ +1 - 51 


u - a ' -0.40 

1 io +2 - 24 
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5 T ( ii ) + p o wer- law 


n 4n+ u a8 
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0.70t[Ji y 6 
0.671^ 
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large error bars. Figure 6 shows the derived radial profiles 
of O/Fe, Ne/Fe, and Mg/Fe ratios. On the whole, the de- 
rived abundance ratios of the O/Fe, Ne/Fe, and Mg/Fe 
are consistent with those within 4 r e , although these ra- 
tios have large error bars. At Suzaku's angular resolution, 
about 30% of the photons from the luminous central re- 
gion within 2 r e escape into the ring of 2-4 r e . As a result, 
the derived abundance ratios of the two regions are consis- 
tent with each other. At 3-6 r e , the signal-to- noise ratio 
is not high enough to constrain the abundance pattern. 

5. Discussion 

5.1. Comparison with abundance patterns in the ISM of 
elliptical galaxies 

Figure 7 shows the pattern of abundance ratios in the 
ISM of NGC 4382. Here, we use those derived from spec- 
tral fitting with the IT model (Table 1). The abundance 
ratios of O/Fe, Ne/Fe, and Mg/Fe have subsolar values 
of ^0.2-0.9 in solar units. O/Fe abundance ratio is obvi- 
ously smaller than those of Ne/Fe and Mg/Fe, although 
the values have large error bars. 

Figure 7 also compares the derived abundance pat- 
tern of NGC 4382 with those of other elliptical galax- 
ies, NGC 720 (Tawara et al. 2008), NGC 1399 and NGC 




Fig. 6. Radial profiles of the abundance ratios of O, Ne, and 
Mg relative to the Fe abundance. 
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Table 3. Radial results of NGC 4382 spectral fittings 



ring 


model 


kT 
(kcV) 


O 

(solar) 


Ne 
(solar) 


Mg 
(solar) 


Fe 
(solar) 


xVd.o.f. 


0-2 r e 


lT+power-law 


n Qi 8+0.030 
u.oio_ 026 


> 0.17 


>0.39 


> 0.17 


> 0.60 


638 / 532 


2-4 r e 


lT+power-law 


n oi 9 +0.024 
u.oiz_ 033 


>0.15 


> 0.30 


>0.41 


>0.49 


631 / 538 


3-6 r e 


lT+power-law 


OO7+0.044 
U.OZ 1 _q.039 


u-^U-o.ig 


> 0.12 


n co+5.45 
U - Oz -0.38 


u - 40 -0.26 


752 / 619 



ring 


model 


O/Fe 
(solar) 


Ne/Fe 
(solar) 


Mg/Fc 
(solar) 


0-2 r e 

2- 4 r e 

3- 6 r e 


lT+power-law 
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1404 (Matsushita et al. 2007b), and NGC 4636 (Hayashi 
et al. 2009). Here, the abundance values in solar units 
are converted using the new solar abundance table of 
Lodders (2003), since the other papers use different so- 
lar abundance tables, those of Feldman (1992) or Anders 
& Grevesse (1989). O/Fe and Ne/Fe abundance ratios in- 
creased by 60% than those adopting the solar abundances 
of Feldman (1992), while the Mg/Fe abundance ratio de- 
creased by 30%, compared with that adopting the solar 
abundances of Anders & Grevesse (1989). 

The abundance ratios of the four elliptical galaxies have 
a similar pattern with small scatter, except for the Ne/Fe 
ratios, and are not very different from the solar abun- 
dances of Lodders (2003). The average ratios of O/Fe, 
Ne/Fe, and Mg/Fe for these four elliptical galaxies are 
0.77, 1.10, and 0.71 in solar units, respectively. As shown 
in Figure 7, the O/Fc ratio of NGC 4382 is smaller than 
those of the elliptical galaxies. Adopting the 99% confi- 
dence error, the upper O/Fe value of NGC 4382 becomes 
0.61 in solar units, which is still smaller than the values of 
the ellipticals. The Ne/Fe ratio of the galaxy is consistent 
with the smallest value for the ellipticals. The uncertainty 
in the spectral modeling of the Fe-L lines may cause the 
scatter in the Ne/Fe ratio in these galaxies and also in the 
ICM in groups and clusters of galaxies (e.g. Matsushita 
2001, Komiyama et al. 2009). The Mg/Fe ratio of the SO 
galaxy has a large error and is consistent with those of 
the ellipticals. Since O is much more abundant in number 
relative to H than Ne and Mg, hereafter we will discuss 
possible differences in the metal enrichment histories of 
the ISM in NGC 4382 and the elliptical galaxies. 

5.2. Contributions from SNe la and SNe II 

The theoretical SNe la and SNe II yields are also plotted 
in Figure 7. Here, we refer to the SNe la yields of the W7, 
WDD1, and WDD3 models calculated by Iwamoto et al. 
(1999). The SNe II yields were taken from Nomoto et 
al. (2006), which is an average over the Salpeter initial 
mass function of stellar masses from 10 to 50 M Q with a 
progenitor metallicity of Z = 0.02. 

The ISM abundance ratios of NGC 4382 and the ellipti- 
cal galaxies are located between those of SNe la and SNe 
II. This means that the metals in the ISM are a mixture 
of SNe la and SNe II yields. 



Since O is predominantly synthesized by SNe II, the 
observed lower O/Fc ratio in the ISM of NGC 4382 indi- 
cates larger contributions from SNe la to the ISM in this 
SO galaxy than in elliptical galaxies. In early-type galax- 
ies, the ISM comes from stellar mass loss and is polluted 
by present SNe la. To explain the low O /Fe ratio in NGC 
4382, the present SN la rate may differ between NGC 4382 
and elliptical galaxies, or stars in this galaxy may contain 
more SNe la products than those in elliptical galaxies. 

5.3. SN la rate of SO and elliptical galaxies 

Assuming that the Fe abundance in the ISM is mainly 
the sum of stellar metallicity and the contribution from 
SNe la, it can be expressed as 

(tLHt)M^)^ 

= (g) +35 ii.°xio-»/yr/L B n°-^; (solar)(1) 

(see Matsushita et al. 2003 for details). Here, (Fc/H)„ 
is the stellar Fe abundance synthesized by SNe la, 6*sn 
is the SN la rate, Mg^ is the mass of Fe synthesized 
by one SN la, a* is the mass loss rate of stars, and 
z soiar * s the Fe metallicity of gas with the solar abun- 
dance. Assuming the age to be 13 Gyr, we used the 
mass loss rate from Ciotti et al. (1991), which is approx- 
imated by 1.5 x 10~ n 3 M@/yi. Here ti$ is the age 
in units of 15 Gyr, and Lb is the B-band luminosity. 
MFe produced by one SN la explosion is likely to be ~ 
0.6M Q (Iwamoto et al. 1999). Using the SN la rate in 
early-type galaxies, which was estimated to be ~0.1-0.3 
SN la/100 yr/10 10 L Q in previous studies (Crane et al. 
1977; Barbon 1978; Cappellaro et al. 1997; Cappellaro et 
al. 1999; Turatto et al. 1999; Sharon et al. 2007; Mannucci 
et al. 2008), the present contribution of SNe la to the Fe 
abundance in early-type galaxies becomes 2.9-8.7 solar. 

Recently, Mannucci et al. (2008) reported that the clus- 
ter membership and morphology of early-type galaxies 
may affect the SN rate. In clusters, the derived SN la rate 
of SO galaxies, 0.121±%%% SN la/100 yr/10 10 Af Q , tends 

to be higher than the O.053±ol81o SN Ia/lOOyr/lO lo Af 
of elliptical galaxies. Since the confidence level of these 
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Fig. 7. Abundance ratios of O, Ne, and Mg against Fe in the ISM of NGC 4382 from the IT model fit (filled red circles). For 
comparison, we plotted the abundance ratios of elliptical galaxies NGC 720 (black open circles; Tawara et al. 2008), NGC 1399 
(blue open squares; Matsushita et al. 2007b), NGC 1404 (purple open triangles; Matsushita et al. 2007b), and NGC 4636 (magenta 
crosses; Hayashi et al. 2009). Solid line represents the abundance pattern of SNc II (Nomoto et al. 2006), and dot-dashed, dotted, 
and dashed lines represent those of WDD1, WDD3, and W7 of SNe la, respectively (Iwamoto et al. 1999). 



values is 1 a, the difference is rather marginal. Assuming 
the stellar B-band mass-to-light ratio of NGC 4382 and 
the elliptical galaxies to be 3.3 and 5.0 Mq/Lq (Nagino 
& Matsushita 2009), respectively, the contributions to Fc 
abundance from SNe la in NGC 4382 and the elliptical 
galaxies are 5.9-20.9 and 4.8-11.9 solar, respectively. As 
a result of the large error bars, these values are consistent 
with each other. 

The Fe abundances expected from the observed SN la 
rate are higher than the Fe abundances in the ISM in 
the early-type galaxies observed with Suzaku. This result 
indicates a low SN la rate, although the absolute values 
of Fe abundances have rather large errors owing to the 
elongated shape of the confidence contours. We note that 
a higher SN rate might be related to younger stellar age, in 
which case the stellar mass loss rate might also be higher. 
For example, assuming a stellar age of 10 Gyr, the mass 
loss rate increases by 1.4 times, and the contribution of 
SNe la to the Fe abundance decreases by 0.7 times. 

5.4- Stellar O/Fe abundance ratio and evolution of SO 
galaxies 

The other possibility is that the lower O/Fe ratio of 
NGC 4382 may be because of a difference in the SN la 
products included in stars in this and the elliptical galax- 
ies. Since a longer star-formation time scale yields more 
SN la products in stars, the difference in O/Fc ratio in 
stars can constrain the star-formation history. Suzaku 



observes the ISM of the entire galaxy and constrains the 
abundance pattern in stars of the entire galaxy, while op- 
tical observations are limited to the very center. 

Optical observations found higher abundances of Fe 
than of a-elements in stars in the central region (within 
2 arcsec) of NGC 4382; [a/Fe] = 0.12 ±0.06 at the 1 a 
confidence level (McDermid et al. 2006). Generally, gi- 
ant elliptical galaxies have higher stellar abundances and 
higher a/Fe ratios. For example, the [a/Fe] of NGC 720 
is 0.37 ±0.05 within r e /8 from the center of the galaxy 
(Humphrey et al. 2006). Therefore, the low O/Fe ratio in 
the ISM of NGC 4382 may reflect a lower [a/Fe] in stars. 
XMM-Newton and Chandra observations show that about 
half of the photons from NGC 4832 within 8 r e are emitted 
from the region within 2 r e (Nagino & Matsushita 2009). 
Therefore, we can conclude that at least within a few r e 
the ISM shows the low O/Fe ratio. This region is much 
larger than that within 2 arcsec, to which the optical ob- 
servations (McDermid et al. 2006) are limited. In other 
words, the Suzaku observations provide evidence that not 
only the stars in the central region of NGC 4382, but also 
those within at least a few r e have a smaller a/Fe ratio 
than those of the four elliptical galaxies. 

A number of authors found a clear correlation between 
the a/Fe ratio and system mass (e.g., Worthcy et al. 
1992; Nelan et al. 2005; Thomas et al. 2005; Graves et 
al. 2007). Thomas et al. (2005) found that galaxies with 
stellar velocity dispersion a > 200 km s _1 have a rela- 
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tively higher [a/Fe] (> 0.2) and also have higher metal- 
licitics. NGC 4382 has a lower ISM temperature of ~0.3 
kcV, and its stellar velocity dispersion a is 179±6 km s _1 
(McDermid et al. 2006). Therefore, the lower O/Fe ra- 
tio in NGC 4382 may reflect a correlation with the sys- 
tem mass. Usually, these metallicity and a/Fe correla- 
tions with system mass are interpreted as showing that 
higher-mass systems have experienced more rapid star 
formation than lower-mass galaxies (e.g. Thomas ct al. 
2005; McDermid et al. 2006). Thus, most stars in larger 
systems would have formed before SNe la began to syn- 
thesize large amounts of Fe. 

Considering the fractional evolution of SO and spiral 
galaxies in clusters, spiral galaxies might have changed 
into SO galaxies (e.g. Dressier et al. 1997; Fasano et al. 
2000; Treu et al. 2003; Kodama et al. 2004; Postman et al. 
2005; Smith et al. 2005; Desai et al. 2007; Poggianti et al. 
2009). If so, stars in SO galaxies would contain more SN 
la products than those in elliptical galaxies, reflecting the 
fact that stars in spiral galaxies have lower [a/Fe] ratios 
than those in giant ellipticals. For example, the value of 
[a/Fe] of the Sun is zero. If SO galaxies in clusters gener- 
ally have lower O/Fe ratios in the ISM, the hypothesis of 
morphology changes would be supported. To investigate 
further, we need more sample galaxies for which to derive 
abundance patterns in the ISM using Suzaku observations. 

6. Summary and conclusion 

With the Suzaku observation of SO galaxy NGC 4382, 
we measured the ISM temperature, metal abundances of 
O, Ne, Mg and Fe, and their abundance ratios for the 
region within 4 r e of the galaxy's center. The temperature, 
0.3 kcV, and O/Fe ratio, 0.3 in solar units, in the ISM of 
this galaxy are smaller than those in elliptical galaxies 
observed with Suzaku, NGC 720 (Tawara et al. 2008), 
NGC 1399 and NGC 1404 (Matsushita et al. 2007b), and 
NGC 4636 (Hayashi et al. 2009). The lower O/Fe ratio of 
the ISM may reflect a higher rate of present SN la or a 
lower a/Fe ratio in stars in NGC 4382 (McDermid et al. 
2006) than in those in the four elliptical galaxies. 

We thank to the referee for a careful reviewing and 
many helpful comments. We also thank all members of 
the Suzaku hardware and software teams and the science 
working group. 
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